Preprint typeset in JHEP style - HYPER VERSION FSU-HEP-020530 



> 



o 

Oh 



Radiative Neutralino Decay in Supersymmetric Models 



(N 
O 

o 

Howard Baer and Tadas Krupovnickas 

m : 

Department of Physics, Florida State University 
Q\ I Tallahassee, FL, USA 32306 

E-mail: 



baerOhep .fsu.edu, tadasOhep . f su . edu 



Abstract: The radiative decay Z 2 — ► Z17 proceeds at the one- loop level in the MSSM. It 



can be the dominant decay mode for the second lightest neutralino Z 2 in certain regions 
of parameter space of supersymmetric models, where either a dynamical and/or kinematic 
(Nl ■ enhancement of the branching fraction occurs. We perform an updated numerical study 

of this decay mode in both the minimal supergravity model (mSUGRA) and in the more 
general MSSM framework. In mSUGRA, the largest rates are found in the "focus point" 
region, where the [i parameter becomes small, and the lightest neutralinos become higgsino- 



like; in this case, radiative branching fraction can reach the 1% level. Our MSSM analysis 
includes a scan over independent positive and negative gaugino masses. We show branching 
fractions can reach the 10-100% level even for large values of the parameter tan (3. These 
^ ■ regions of parameter space are realized in supergravity models with non-universal gaugino 

masses. Measurement of the radiative neutralino branching fraction may help pin down 
underlying parameters of the fundamental supersymmetric model. 



Keywords: Bupersymmetry Phenomenology, Supersymmetric Standard Model, Rare 
Decays. 



1. Introduction 

One of the major goals of collider experiments is to verify or disprove the existence of 
weak scale super symmetric matter [Q]. Given sufficient energy and integrated luminosity, 
sparticles should be produced at large rates, if indeed they exist. Once produced, sparticles 
will decay via a cascade into other sparticles until the state containing the lightest SUSY 
particle (LSP) is reachedQ. If kinematically accessible, tree level two-body decays will 
comprise the dominant sparticle branching fractions. If two-body modes are suppressed or 
closed, then tree level three-body decays can dominate. It is also possible for sparticles to 
decay into two-body final states where the decay is mediated by loop diagrams. Important 
examples include g — > gZi, Z{ — > Zj-f (i, j = 1 — 4 with i > j) and t\ — > cZ\ for sparticles, 
and h — > 77 for the lightest SUSY Higgs scalar. 

The radiative decay Z\ — * Zjj has been examined previously in a number of papers, 
since it can give rise to a distinctive energetic, isolated photon plus missing energy signature. 
The decay proceeds through loops containing intermediate charged scalars and charginos. 
After several approximate calculations^], the complete calculation was presented in both 
the non-linear and linear gauge by Haber and WylerQ. These authors note that the 
radiative neutralino decay branching fraction may be large in MSSM parameter space 
regions where M\, M2, \fi\ <C Mz (light neutralino radiative decay), or when \/j,\,Mz "C 
Mi, M2 ~ TeV (higgsino to higgsino radiative decay). 

Detailed numerical studies of the Z2 — > Z17 decay rate have been performed by Am- 
brosanio and Mele||]. These authors found that the Z2 — > Z17 decay rate can in fact be 
the dominant Z2 decay mode in regions of MSSM model parameter space where tree level 
two-body Z2 decay modes are not allowed, and the competing Z2 three-body decay modes 
are dynamically and/or kinematically suppressed. 

At low tan/3, the three-body neutralino decay Z2 — ► Z\j f occurs dominantly via 
five Feynman diagrams involving virtual sfermion and Z boson exchange. The dynamical 
suppression of the three-body decay modes occurs when one of Z2 and Z\ is gaugino-like, 
while the other is higgsino-like. Since the Ziff vertex contains the gaugino component of 
the neutralinos, the decay via a virtual sfermion /* is suppressed if one (or both) of Z2 or Z\ 
is higgsino-like. In addition, the ZZ2Z\ coupling contains the higgsino components of the 
neutralinos so that if one of Z\ and/or Z2 is gaugino-like, then decay via the Z* diagram is 
suppressed. The region of MSSM parameter space where dynamical suppression of three- 
body decay modes occurs is when the gaugino masses M\ ~ M2, tan/? ~ 1 and [i < with 
\fi\ being sufficiently small. In light of recent limits from LEP2 that > 103.5 GeV and 
nih > 114.1 GeV (for a SM-like h boson), the region of dynamical suppression is largely 
excluded, since mh is quite light at very low values of tan/?. 

In addition, there can exist a region of reduced dynamical suppression. In this case, 
three-body decays via Z* can be suppressed if one (or both) of the neutralinos is gaugino- 
like (thus suppressing Z2 decay via Z*), while the sfermions are quite heavy (thus sup- 
pressing decay through /*). This interesting case occurs naturally as we shall see in the 
"focus-point" scenario^] of the minimal supergravity (mSUGRA or CMSSM) model[^|. 

Finally, there can exist regions of kinematical suppression of three-body Z2 decays. 
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In this case, mj? 2 ~ rug so that Z2 — ► Z\f f is suppressed kinematically by a factor e 5 , 
where e = (1 — /m^ ). The radiative decay Z2 — > Z17 is suppressed by e 3 , so that 
as the m^ 2 — mass gap approaches zero, the radiative decay increasingly dominates. 
The kinematic suppression region occurs when tan (3 ~ 1 or when M\ ~ M2 . In the region 
of parameter space where kinematical suppression occurs, however, the small mass gap 
between neutralinos also means that the final state photon energy will be small, unless 
enhanced by a Lorentz boost if Z2 is produced either directly or via cascade decay with a 
high velocity. 

The radiative neutralino decays have received some recent attention by Kane et al.|| 
in attempts to explain an ee77 event seen by the CDF experiment in run 1 of the Fermilab 
Tevatron|J. These "higgsino- world" scenarios seem to need tan/3 ~ 1, /i < and a light 
top squark. These scenarios seem largely excluded now by LEP2 limits on the masses of 
the Higgs bosons, which require tan j3 > 2 — 3. 

In this paper, we re-examine the radiative neutralino decay branching fraction as a 
function of supersymmetric model parameter space. Our calculations include a number of 
improvements over previous analyses. 

1. We implement the complete Z2 — * Z17 decay calculation into the event genera- 
tor IS A JET v7.64[10]. Our calculation includes full third generation mixing effects 
amongst the squarks and sleptons. 

2. Our treatment of neutralino three-body decays includes the effect of all third gener- 
ation Yukawa couplings and sfermion mixings in the decay calculations^]. These 
are especially important at large values of the parameter tan (3, where the Yukawa 
couplings fb and f T can become large. In addition, along with the five Z2 — ► Z\f f 
decay diagrams proceeding via sfermion and Z boson exchange, we include decays 
through h, H and A bosons. These diagrams can be important, especially at large 



We include the latest constraints from LEP2. As noted before, these constraints 
eliminate many of the regions of parameter space where Z2 — ► Z17 can be large, 
especially around tan/3 ~ 1. 

We present results for the paradigm mSUGRA model. We also present a scan of 
MSSM parameter space which includes negative as well as positive gaugino masses. 
We find the radiative decay branching fraction can be large in a particular SUGRA 
model with non-univeral gaugino masses. 

Finally, in our numerical scans, we take note of parameter space regions where the 
radiative neutralino branching fraction is not dominant, but may be nonetheless 
observable. For instance, at the CERN LHC, many millions of sparticle production 
events may be recorded each year, provided sparticles are light enough. In this case, 
branching fractions of a few per cent or even less may be interesting, and allow one 
to measure the neutralino radiative decay branching fraction. Since the Z2 — * Z-\pj 
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branching fraction is very sensitive to parameter space, its measurement may help 
determine or constrain some of the fundamental SUSY model parameters. 



The rest of this paper is organized as follows. In Sec. |2], we examine the radiative neu- 
tralino branching fraction in the mSUGRA model. We find it can reach levels approaching 
1%, especially in the well-motivated focus-point region. In Sec. |3|, we examine the radiative 
decay branching fraction in the MSSM. In this case, we find regions of parameter space 
where the branching fraction can approach 50-100%, even for large values of tan/3. These 
regions occur dominantly where M\ ~ M 2 or when M\ ~ —M 2 . We note in Sec. |I| that 
the latter MSSM regions of enhancement naturally occur in a particular SUGRA model 
with non-universal gaugino masses. We map out the interesting regions of parameter space 
where the radiative decay branching fraction can be large. Finally, in an appendix, we 
present our formulae for radiative neutralino production, cast in a form consistent with the 
ISAJET event generator program. 

2. Z 2 — > Zi'j branching fraction in mSUGRA model 

Our first goal is to examine the radiative neutralino branching fraction in the paradigm 
mSUGRA (or CMSSM) model. Most analyses of supersymmetric phenomena take place 
within this model, so it is convenient to know the rates for radiative neutralino decay in 
relation to other expected phenomena. In the mSUGRA model, the MSSM is assumed 
to be valid between the grand unified scale (Mqut ~2x 10 16 GeV) and the weak scale 
(~ 1 TeV). At Mqut, it is assumed all scalar masses are unified to the parameter mo, 
while all gaugino masses unify to m\i 2 - In addition, trilinear soft SUSY breaking masses 
unify to Aq. The gauge and Yukawa couplings and soft SUSY breaking parameters evolve 
from Mqut to M wea ^ according to the well-known MSSM renormalization group equations 
(RGEs). At M wea k it is assumed electroweak symmetry is broken radiatively, which fixes 
the magnitude of the superpotential [i parameter, and allows one to effectively trade the 
bilinear soft breaking parameter B for the ratio of Higgs vevs tan (3. Thus the parameter 
space is characterized by 

m 0> m i/2> Aq, tan/3 and sign(fi), (2-1) 

where we take the top quark pole mass mt = 175 GeV. In our mSUGRA model calculations, 
we utilize the iterative RGE solution embedded within the program ISAJET v7.64[10]. This 
version also includes our radiative neutralino decay calculation. 

In Fig. [[], we show the mo vs. rrii/ 2 plane for tan/3 = 10, Aq = and /i > 0. The 
red-shaded regions are excluded either by the presence of a stau LSP (upper left) or a lack 
of REWSB (lower and right). The pink shaded region is excluded by LEP2 searches for 
charginos and a light Higgs scalar. The upper and left regions indicate where Z 2 two-body 
decays occur: either to l\i (£ denotes any of the charged or neutral leptons) on the left, or 
to Z\Z or Z\h in the upper plane. The remaining regions allowed by LEP2 have dominant 
three-body Z 2 decays. The green shaded region has 0.001 < BF(Z 2 — ► Z17) < 0.005 while 
the yellow has 0.005 < BF(Z2 — > Z17) < 0.01. In these regions, where ~ 250 GeV, 
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we also have rrig ~ 650 GeV. Then gg cross sections exceed ~ 2000 fb at the CERN LHC, 
leading to at least 20,000 SUSY events per 10 fb -1 of integrated luminosity. After suitable 
cuts are made on number of jets, jet Er and Et, a nearly pure sample of SUSY events 



will remain] 14 1 . These SUSY candidates can then be scrutinized for the presence of hard, 
isolated photons, and may allow a determination at some level of the neutralino radiative 
decay branching fraction. 

We note that the regions of mSUGRA parameter space with the largest Z2 — > Z-\pj 
branching fraction occur at rather large mo and low values. The very largest allowed 
mo areas correspond to the focus-point region of Feng, Matchev and MoroiQ, where fine- 
tuning may be low, but squark and slepton masses are at the multi-TeV level. These 
heavy sfermion masses may be sufficient to suppress many flavor-changing or CP violating 
processes that could arise due to the presence of off-diagonal flavor changing soft SUSY 
masses, or imaginary components of soft SUSY breaking terms. 

We also show the regions of Z2 — ► Z\"f branching fraction for large tan f3 = 45 in Fig. |2[ 
Again, we see a significant part of the three-body Z2 decay region inhabited with possibly 
measurable branching fractions for radiative neutralino decay. The region with the largest 
branching fraction is adjacent to the right-side region forbidden by lack of REWSB. In 
that region, the fi parameter is becoming small, and the light neutralinos are becoming 
increasingly higgsino-like. In addition, sfermion masses are heavy, so we are in a region 
of reduced dynamical suppression of three-body decay modes, also with some kinematic 
suppression. 

In both the low and the high tan/3 cases, it is important to note that the Z2 — > Z\^ 
branching fraction reaches no greater than 1%, so is never dominant. However, regions do 
occur where the radiative decay should be detectable, and the branching fraction should 
be measureable. Detailed simulations including detector effects for identifying isolated 
photons will be necessary to ascertain the precision with which such measurements can be 
made. 

3. Z 2 — > Zi7 branching fraction in the MSSM 

It is apparent from the numerical studies of Ambrosanio and Mele[H] that the greater 
freedom of parameter choices in the unconstrained MSSM can allow for regions of parameter 
space with larger Z2 — > Z\j branching fractions than in the mSUGRA model. In this 
section, we present a scan over unconstrained MSSM parameters. To make the analysis 
tractable, we do assume all MSSM soft SUSY breaking parameters and Yukawa couplings to 
be purely real, and we take all sfermion masses to be degenerate, with no intergenerational 
mixing. For simplicity, we take At = Af, = A T , and fix nig = 1500 GeV. Our scan of MSSM 
parameter space is then over the following region: 

-1500 GeV < Mi < 1500 GeV, 

< M 2 < 1500 GeV, 

-1500 GeV < n < 1500 GeV, 

-2000 GeV < A t < 2000 GeV, 
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D.DD5< BF<D.D1 « 0.001< BF <D.DD5 
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Figure 1: Regions of mSUGRA model parameter space with significant Z 2 — > Z17 branching 
fraction, in the m us. mi/2 plane for tan/3 = 10, A§ = and /i > 0. The red shaded regions are 
theoretically excluded while pink regions are excluded by sparticle searches at LEP2. The green 
region has 0.001 < BF(Z 2 -> Zij) < 0.005, while the yellow region has 0.005 < BF(Z 2 -> Zfy) < 
0.01. 

150 GeV < m A < 1500 GeV, 

2 < tan/3 < 50, 
100 GeV < m f ~ < 1500 GeV, 

where we have the freedom to choose one of the gaugino masses, in this case M2, to be 
real and positive. We require mjj > 103.5 GeV and > 114.1 GeV, according to LEP2. 
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D.DD5< BF<D.D1 • D.001< BF<D. DOS 
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Figure 2: Same as Fig. 1, except tan/3 = 45. 

The latter Higgs mass requirement is overly harsh in that the limit on the lightest SUSY 
Higgs boson is somewhat lower, but this doesn't alter our conclusions. We also require a 
mass gap — > 5 GeV, so that the photon that appears in scattering events from 
any of these models has a detectable energy. 

Our results are presented in Figs. || and ||, as plots of branching fraction {BF) versus 
model parameter. Fig. ||a and show the branching fractions versus the gaugino masses 
Mi and M%. We see immediately that in fact there do exist models with substantial 
radiative decay branching fractions in the range 10%-100%. For both gaugino masses, these 
large branching fraction models populate the region with \Mi\, M<i < 1000 GeV, with the 
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Figure 3: Plot of BF(Z 2 -> Z17) versus the parameter a) Ml, 6) M 2 , c) ^ and cQ A t in the MSSM 
model. 

distribution of points being somewhat asymmetric for M\ being positive or negative. If M\ 
becomes large negative or M2 becomes large, the branching fractions can remain large, but 
the mass gap m^ 2 — decreases below 5 GeV. However, for large positive M\ the mass 
gap stays large, but in this case we do not get large branching fractions because either 
two-body decay modes of Z2 turn on, or else Z2 — ► W\f f becomes large, thus supressing 
the radiative decay branching fraction. The models with \Mi\, M2 near zero are usually 
excluded by constraints from LEP onZ-> ZiZj decay or from LEP2 on the chargino mass. 

Fig. ||c shows the BF versus the [i parameter. Again, small values of are excluded 
by LEP and LEP2 constraints. However, we see that really large radiative decay branching 
fractions only occur for very large values of ~ 300 — 400 GeV. The suppresion of lower 
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Figure 4: Plot of BF{Z,2 — > ^17) versus the parameter a) itla, b) tan/3, c) mi and ef) mg 2 — m z l 
in the MSSM model. 

values of comes from the mass gap requirement. Fig. shows the branching fraction 
versus the parameter A t . Here, large values of branching fraction are possible for all values 
of At, although the very largest seem to favor large values of \At\. 

In Fig. H| we show the BF versus the parameters uia, tan (3 and mj in frames a, b and 
c. We see that large branching fractions are possible for all values of pseudoscalar mass 
rriA, and for all values of tan {3. Large values of BF have some preference for mj ~ 300 — 400 
GeV, so that two-body decays Z2 — > // are closed, and three-body decays Z2 — > Z\f f 
through an intermediate / have some suppression. 

In the last frame, Fig. we show the branching fraction not against an input 
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parameter, but against the neutralino mass gap, — . Here it is clear that the 
largest values of BF prefer the smallest mass gaps, where there exists at least some amount 
of kinematical suppression of three-body Z2 decays. The large branching fractions have a 
cut-off around 90-100 GeV, where Z 2 — > Z\Z decays become possible. 

To gain more understanding of the parameter space where radiative decays are large, 
we plot in Fig. || only models with BF > 10% (black dots) or > 1% (red dots), in 
correlated parameter planes. Again, we require m^ 2 — > 5 GeV, and all LEP and 
LEP2 constraints to be satisfied. From Fig. |5]a in the M\ vs. M 2 plane, we see that models 
with the largest BF require \M\\ ~ M 2 , a feature noted by Ambrosanio and Mele in their 
analysis ||. The region with M\ ~ M2 turns out to be quite narrow, while the region with 
Mi ~ —M2 has slightly more breadth to it. This latter point is even more apparent when 
we relax to the condition BF > 1%. 

In Fig. |5]& and [|c, we show models in the M\ vs. fi and M2 vs. \x planes. For models 
with BF > 10%, evidently \M\\ < and also M2 < If we relax to models with 
BF > 1%, then these requirements no longer hold. In frame [5]d, we show models in the 
M 2 /Mi vs. [i plane. Models with BF > 10% clearly favor M 2 /Mi ~ 1, with either sign 
of [i being equally preferable. If we allow models with BF > 1%, then the M\ ~ ±M 2 
restriction can be relaxed, and rather large ratios of M2/M1 are allowed, especially if 
is small. 

An important element of our analysis is the size of the Z2 — Z\ mass gap, so there is 
sufficient photon energy for observations. In Fig. ||, we show the mass gap m^ 2 ~ m z 1 
versus the parameters Mi, M2 and tan/3 in frames a-d, respectively. From Fig. ^a, we 
see that models with BF > 10% prefer large absolute values of the \i parameter, but for 
these models the mass gap is bounded typically by 40-60 GeV. By relaxing to models with 
BF > 1%, we pick up models with mass gaps ranging to 100 GeV. However, for models 
with small values of where higgsino^ higgsino transitions are possible, the mass gap 
is very small (~ 20 GeV), and these cases will be more difficult to observe. Figs. ^6 and 
|6]c show that the mass gap for models with BF > 10% is typically less than 30-40 GeV 
for \Mi\, M 2 ranging up to 1000 GeV. The models with BF > 1% and |Mi|, M 2 > 1000 
GeV typically have a small mass gap, less than 10-15 GeV. Finally, in Fig. [6|d, the mass 
gap versus tan/3 shows no major correlations, and significant mass gaps can be achieved 
for any tan (3 values. 

To conclude this section, we summarize our findings by noting that the largest neu- 
tralino radiative decay branching fractions are found in regions of parameter space with 

I Mi I, M 2 ~ 100 - 1000 GeV, 
\pi\ > \Mi\, M 2 and 
|Mi| ~ M 2 . 

4. Z2 — > Zi'-f branching fraction in a SUGRA model with non-universal 
gaugino masses 

So far, our findings are that in the mSUGRA model the radiative decay BF never exceeds 
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Figure 5: Plot of MSSM model points with BF(Z 2 -> Z17) > 0.1 and - > 5 GeV in a) 
Mi vs. Mi space, b) M\ vs. fi space, c) M2 vs. /i space and d) M2/M1 vs. [i space. 



1%, and is largest in regions with large mo such as the focus point region. For the MSSM 
case, we found distinct regions of parameter space where the radiative branching fraction 
can range up to 10%-100%, although these conditions are never realized in the mSUGRA 
model. In this section, we note that there do exist models where the above MSSM param- 
eter space restrictions are realized, namely in SUGRA models with non-universal gaugino 
masses. 

Since supergravity is not a renormalizable theory, in general we may expect a non- 
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Figure 6: Plot of the mass gap to^ 2 — versus a) /i, b) Mi and c) M2 for MSSM model points 
with TSF(Zi — > Z17) > 0.1 and mg 2 — wi^ > 5 GeV. In rf), we plot the mass gap mg 2 — versus 
the parameter tan (3 in the MSSM. 



trivial gauge kinetic function. Expanding the gauge kinetic function as 

Iab = Sab + h AB /M P + ..., (4.1) 

where the fields Hab transform as left handed chiral superfields under supersymmetry 
transformations, and as the symmetric product of two adjoints under gauge symmetries, 
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Table 1: Relative gaugino masses at Mqut and Mz in the four possible irreducible represen- 
tations. 

we parametrize the lowest order contribution to gaugino masses by, 

where the Wa are curl superfields, the are the gaugino fields, and F^ab are the auxiliary- 
field components of the Kab that acquire a SUSY breaking vev. 

If the fields Fh which break supersymmetry are gauge singlets, universal gaugino masses 
result. However, in principle, the chiral superfield which communicates supersymmetry 
breaking to the gaugino fields can lie in any representation in the symmetric product 
of two adjoints, and so can lead to gaugino mass terms that (spontaneously) break the 
underlying gauge symmetry. 

In the context of SU(5) grand unification, F^ belongs to an SU(5) irreducible repre- 
sentation which appears in the symmetric product of two adjoints: 

(24x24) symmetric = 1 24 75 200 , (4.3) 

where only 1 yields universal masses. The relations amongst the various GUT scale 
gaugino masses have been worked out by Anderson |l5fl. The relative GUT scale SU(3), 
SU(2) and U(l) gaugino masses M3, M2 and M\ are listed in Table Q along with the 
approximate masses after RGE evolution to Q ~ Mz- These scenarios represent the 
predictive subset of the more general (and less predictive) case of an arbitrary superposition 
of these representations. The model parameters may be chosen to be, 

mo, M3, Aq, tan/3 and sign(fi), (4-4) 

where M° is the SU(i) gaugino mass at scale Q = M GUT . M% and M° can then be 
calculated in terms of M3 according to Table [l], and the weak scale SUSY particle spectrum 
can again be calculated using ISAJET v7.64. 

We note that the 75 model leads to weak scale gaugino masses in the ratio M2 : 
Mi ~ 6 : —5, which is approaching the conditions required in Sec. || for models with large 
radiative branching fractions. We adopt the GUT scale gaugino mass boundary conditions 
of the 75 model, and explore the model parameter space for large radiative decay branching 
fractions. 

Our first results are shown in Fig. |7| in the M3 vs. mo plane, for tan/3 = 5, Aq = 
and fi > 0. Again, the red shaded regions are excluded theoretically, while the pink regions 
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are excluded by constraints from LEP2. The region with mass gap — > 5 GeV 
is indicated between the dotted contours. The substantial green region has branching 
fractions BF : 1% — 10%, while the yellow region has branching fractions in the range 
BF : 10% — 50%. Almost all the large branching fraction region has mass gap m z 2 ~ m z 1 < 
5 GeV. The value of in the yellow region is generally much smaller than \M\\ and M2, so 
that Z2 — > Z17 consists of a higgsino — > higgsino transition, with a rather small mass gap 
between Z2 and Z\\ there is significant kinematic suppression of three-body decay rates. In 
the lower right unshaded regions, the value of [i is comparable or larger than \M\\ and M2, 
and the Z^-Z\ mass gap is usually larger, so three-body decay rates increase, suppressing 
the radiative decay branching fraction. An exception occurs near mo ~ 3000 GeV and 
M3 ~ 300 GeV where the light neutralinos are transitioning between higgsino-like and 
gaugino-like. In this area, the mass gap drops to very small values of order 1 GeV, and the 
radiative decay is momentarily enhanced. 
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Figure 7: Regions of the 75 SUGRA model parameter space with significant Z 2 — > Z17 branching 
fraction, in the mo vs. m.1/2 plane for tan/3 = 5, Ao = and /1 > 0. The red shaded regions are 
theoretically excluded while pink regions are excluded by sparticle searches at LEP2. The green 
region has 0.01 < BF(Z 2 -> Z17) < 0.1, while the yellow region has 0.1 < BF(Z 2 -> Z17) < 0.5. 

In Fig. ||, we show again the 75 model, but this time with tan [5 = 10. As tan f3 is 
increased beyond 5, the allowed parameter space rapidly decreases. In this case, only a 
small region in the center of the plot gives a viable SUSY spectrum satsifying all constraints. 
In this case, the yellow region represents 1% < BF < 5% (lower values than in Fig. 0). 
The larger tan/3 value allows for a larger mass gap between Z2 and Z\ for small \i values 
compared to Fig. [?] ("i^ 2 ~ m z 1 > ^ GeV throughout the tan/3 = 10 plane). Consequently, 
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Figure 8: Same as Fig. 6, except tan/3 = 10, and now the green region has 0.001 < BF(Z,2 — > 
Zxi) < 0.01, while the yellow region has 0.01 < BF(Z 2 -> Z17) < 0.05. 

three-body decays suffer less kinematical suppression, and the radiative branching fraction 
is lower. 
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We have also examined briefly two cases of "more minimal" supersymmetric models [16]. 
In these models, also known as inverted mass hierarchy (IMH) models, there can exist 
non-universal scalar masses, leading to multi-TeV first and second generation scalars (to 
suppress flavor and CP violating SUSY processes), while third generation scalars remain 
in the sub-TeV regime (as required by naturalness constraints) . The two models examined 
were the radiative inverted hierarchy model (RIMH) Jl7| and the GUT scale IMH model 
(GSIMH)(T^]. In both these models, usually two-body Z 2 decays are allowed. In the 
RIMH model, for some parameter regions the \i parameter can become very small. In that 
case, results are similar to the SUGRA model focus point region already presented with 
branching fractions near the ~ 1% level. 



5. Conclusions 



We have re-calculated the decay rate for radiative neutralino decay within the MSSM. We 
have also made detailed scans of parameter space in the mSUGRA model, the MSSM and 
a SUGRA model with non-universal gaugino masses. In the mSUGRA model, radiative 
neutralino branching fractions for Z 2 can reach up to 1%. While this rate comprises only 
a small fraction of the Z 2 decay modes, it may be sufficiently large to be measureable at 
collider experiments at the CERN LHC, where many g and q events can be produced, and 
Z2 will occur in cascade decays. The rate is largest in the focus point region of parameter 
space. If the branching fraction is observed, then it may aid in the determination of 
underlying model parameters. 

In the MSSM, the Z2 — > Z17 decay rate can reach nearly 100%. In this case, the 
gaugino masses \M X \, M 2 ~ 1000 GeV, with |Mi|, M 2 < \n\ and |Mi| ~ M 2 . The 
branching fraction can be large for all values of tan /3. 

The conditions on parameters of the MSSM for a large radiative neutralino branching 
fraction turn out to be realized in a SUGRA model with non-universal gaugino masses, 
where the hidden sector field Hab transforms as a 75 dimensional symmetric product of 
two adjoints of SU(5). At low tan/3 ~ 5, in this case the branching fraction can reach up 
to 50%, although the mass gap — mj^ tends to be rather small - of order a few GeV. 
For higher tan j3, the parameter space in the 75 model rapidly diminishes. However, scans 
made for the tan/3 = 10 case show the radiative branching fraction reaching up to ~ 5%. 

Our final conclusion is that once a significant sample of SUSY particle candidate 
events is discovered at collider experiments, then that sample should be scrutinized for the 
presence of additional isolated photons, beyond the level expected from bremsstrahlung, 
misidentified jets, etc. If the events contain in addition isolated photons, then it may be 
possible to test whether they come from radiative neutralino decay, and to in fact measure 
that particular branching fraction. If the Z 2 — > Zyy branching fraction is in fact mea- 
sured, then its value should aid in the attempt to extract the values of underlying model 
parameters. 
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Appendix 



In this appendix, we present our formulae for the radiative neutralino decay rate. The 
complete radiative decay width of Zj was calculated by Haber and WylerQ. We organize 
our calculation in a similar fashion to theirs, working in the non-linear gauge. We 
use the Lagrangian conventions of Tata|l|, which is useful for inclusion into ISAJET. The 
relevant Feynman diagrams labelled a-p are presented in Ref. H, and will not be repeated 
here. We have checked that our formula reproduces analytically that of Haber and Wyler 
after converting conventions. 

In addition, we have compared our numerical results for various diagram subsets with 
a calculation by Ambrosanio and Mele which is embedded in the SUSYGEN 3.00/36 
program[ IS]. Again, we find excellent agreement, except for the case of sfermion loop 
diagrams, where mixing effects cause our result to differ from those of SUSYGEN. If we 
neglect mixing effects, then we agree with the SUSYGEN result. 

The radiative decay width is given by 



r(Zj - z a ) 



9 Z ] Zg 
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Here 
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with the separate contributions given by the following: 
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Here Cf = 1 for leptons and Cf = 3 for quarks. 
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9 W4 , 
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2 , (i) 
3v^ 4 



(i) 

For the neutralino mixing element Va , a = 1 if / is an up type quark while a = 2 for 
leptons and down type quarks. 

The g and g' are the SU(2) and U (l)y gauge couplings, while the Uj- are neutralino 
mixing elements, as defined by Tata|l]]. The // is the Yukawa coupling for fermion /, and 
Of is the sfermion mixing angle. The integrals /, I 2 , J and K are defined in Haber and 
Wyler§. 

Also, 
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'» i - -/= (0^3 + 9 v\ 'J sin 7/j - cos 7^ 
-y= [gv 3 +9W) cos 7l + gv\ SU17L, 

-7= + 5 u 4 J sm 7L + cos 7L 



Finally, 



9mnop 



(5.6) 



Here (7™ nop (J = 0) are obtained from (7!-^ by making the replacement sin/3 — > 
— cos 0, cos /3 — > sin j3. In the above, 7^ and 7r are chargino mixing angles as defined by 
Tata; in that workQ, the 9 X , 9 y , B% and 6^ are also defined. 
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